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ABSTRACT 

We present an analysis of X-ray, Ultraviolet and optical/near-IR photometric data 
of the transitional millisecond pulsar binary XSS J12270-4859, obtained at different 
epochs after the transition to a rotation-powered radio pulsar state. The observations, 
while confirming the large-amplitude orbital modulation found in previous studies 
after the state change, also reveal an energy dependence of the amplitudes as well as 
variations on time scale of months. The amplitude variations are anti-correlated in the 
X-ray and the UV/optical bands. The average X-ray spectrum is described by a power 
law with T index of 1.07(8) without requiring an additional thermal component. The 
power law index T varies from ~1.2 to ~1.0 between superior and inferior conjunction 
of the neutron star. We interpret the observed X-ray behaviour in terms of synchrotron 
radiation emitted in an extended intrabinary shock, located between the pulsar and 
the donor star, which is eclipsed due to the companion orbital motion. The G5 type 
donor dominates the UV/optical and near-IR emission and is similarly found to be 
heated up to ~ 6500 K as in the disc state. The analysis of optical light curves gives a 
binary inclination 46° < i < 65° and a mass ratio 0.11 < q < 0.26. The donor mass is 
found to be 0.15 < M 2 < 0.36 M 0 for a neutron star mass of 1.4M 0 . The variations 
in the amplitude of the orbital modulation are interpreted in terms of small changes 
in the mass flow rate from the donor star. The spectral energy distribution from radio 
to gamma-rays is composed by multiple contributions that are different from those 
observed during the accretion-powered state. 

Key words: Interactive binaries - Stars: individual: XSS J12270-4859,1FGL J1227.9- 
4852, 2FGL J1227.7-4853, 3FGLJ1227.9-4854 - gamma-rays: stars- X-rays: binaries - 
Accretion 


1 INTRODUCTION 

Millisecond pulsars (MSPs) are believed to be formed in bi¬ 
nary systems containing old neutron stars (NSs), which are 
spun-up to very short periods during a previous Gyr-long 
phase of mass accretion from an evolved companion. During 
this prolonged phase the binary appears as Low-Mass X-ray 
Binary (LMXB) and when mass accretion ceases a pulsar 
powered by the magnetic field rotation turns on and ap¬ 
pears as a rad i o and gamma-ray pulsar. The MSP recycling 
scenario (lAlpar et al Il982l ) found observational support 
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throu gh the discovery of first MS radio pulsar dBacker et all 
11982h and has drammatically been confirmed by the detec¬ 
tion of few hundreds of galactic radio MSPs in binaries. Fur¬ 
thermore a few tens of accreting N Ss spinning at a few ms 
are known (IPatruno fe Watts! 120121 ). which are believed to 
be progenitors of MSP binaries. The transition between the 
two states was first testified by PSR J1023+0038, a (1.67ms) 
MSP found to be in an accr etion disc-state between 2000- 
2001 (lArchibald et al.l 120091 1. This object was considered 
the ” missing-link” between LMXBs and MSPs. It surpris¬ 
ingly turned on again in the X-ray band in 2013, showing 
X-ray pulsations at the NS spin period and disc emission 
lines in the optical band, thus testifying the presence of 


© 2015 RAS 















2 D. de Martino et al. 


accre t ion at some le vel (|StaBi)ers_et_alJ 201^_ Thkata et al.l 


120141 ; IPatruno et al.l 12014 Bogdanov et al. 20151 ). Another 
striking case is the one of IGR J1825-2452, a binary in the 
M28 globular cluster disco vered as an accretin g MSP dur¬ 
ing an outburst in 2013 (IPapitto et al.l 120131). but iden - 
tified as a radio pulsar few years before I Beginl [2006). 
This source, turned into a radio pulsar at the end of 
the outburst, demonstrating that transitions between the 
accretion-powered and rotation-powered states can occur 
on timescales much shorter than secular evolution. The 
hard X-ray source XSS J12270-4859 (henceforth XSS J1227) 
was identified as an unusual low-luminosity (Lo. 3 _iok e v ~ 
6 x 10 33 erg s _1 for a distance of 1.4 kpc) LMXB, surprisingly 
associated to a high en e rgy gamma-ray Fermi - LAT source 
llde Martino et ahl l20ld . 1 20131 ; iHill et ahl l201ll henceforth 
dM10,dM13). The peculiar X-ray and gamma-ray emission 
from this source was explained in terms of a system host¬ 
ing a fast rotating NS in a propeller state llPapitto et al.l 
120141 ). In late 2012/beginning 2 013 it faded to the low- 
est X-rav dBogdanov et ah 20141) and optical dBassa et al.l 
120141 ; Ide Martino et al.l |2014|~ henceforth dM14) levels ever 
observed. Radio foll ow-ups revealed a fast (1.69 ms) ro¬ 
tating radio pulsar (IRov et al.l 120151 ). Pulsed gamma-rays 
were also recently detecte d during the rotation-powered 
state dJohnson et al.l i2015l l. X-ray pulses at the NS spin 
coul d then be detec ted in the previous disc-accretion state 
dPapitto et alj|2015h . but not in the radio pulsar state be¬ 
cause it is too faint. These three examples open the chal¬ 
lenging possibility that during the evolution there is a phase 
where the NS ’’swings” back and forth from a radio pul¬ 
sar to a sub-luminous (Lx < 10 34 ergs _1 ) disc accretion 
state. Analysing the spin period distribution of the accret¬ 
ing LMXBs, th e rad io MSPs and those undergoing Type-I 
bursts, IPapitto et aid (120131 ) argued that nuclear MSPs are 
in an earlier evolutionary stage, while accreting MSPs are 
in a phase much closer to the switch-off and the turn-on of 
a radio MSP. Thus, the transitions between the two states 
may reflect the interplay between the NS magnetosp here and 
the mass transfer rate from the donor star (see Stella et aid 
Il994l ; ICampana et aid 1 19981 : iBurderi et aid 12001 ). It is also 
argued that some, possibly most, of the so-called ’’redback” 
MSP binaries, which have non-degenerate dono rs, can dis¬ 
play such state transitions dRoberts et aid [20 15). Thus, all 
redbacks are potential ’’swinging” MSP s and searches for 
analo gous sources are ongoing (e.g. iBogdanov fc Halpernl 
120151) . 


The two known galactic field redbacks, PSR J1023+0038 
and XSS J1227, linger in a sub-luminous disc accreting state 
since mid-2013 and in a radio pulsar state since end-2012, 
respectively. They represent ideal study cases to understand 
how transitions occur in redbacks. We here present the anal¬ 
ysis of observations in the X-rays and UV band of XSS J1227 
during the current rotation-powered state acquired with 
XMM-Newton in 2014 June and in the optical with the REM 
telescope in 2015 January and February aiming at studying 
the 6.91 h orbital variability of the source. We also compare 
the X-ray and UV/optical variability observed in 2013 De¬ 
cember whe n it was first found to h ave transited to a disc- 
free state bv IBogdanov et ahl d2014l ). lBassa et al . Idml) and 
dM14. We further analyse the spectral energy distribution 
(SED) from radio to gamma-rays to identify the emission 
components in the radio pulsar state. 


Table 1. Summary of the observations 


Telescope 

Instrument 

UT Date 

UT Start 

T 

expo 

(ks) 

XMM-Newton 

EPIC-MOSl 

2014-06-27 

00:54:49 

41.8 


EPIC-MOS2 

2014-06-27 

00:55:23 

41.8 


OM-U 

2014-06-27 

01:03:50 

42.5 


Telescope 

Instrument 

UT Date 

UT Start 

T 

expo 

(sec) 

# 

REM 

ROSS2 

2015-01-20 

02:46 

300 

30 



2015-01-21 

02:40 

300 

31 



2015-01-23 

02:32 

300 

37 



2015-02-12 

01:16 

300 

43 



2015-02-13 

01:24 

300 

44 



2015-02-14 

01:10 

300 

51 


REMIR 

2015-01-20 

02:40 

300* 

32 



2015-01-21 

02:40 

300 

31 



2015-01-23 

02:24 

300 

38 



2015-02-12 

01:20 

300 

44 



2015-02-13 

01:10 

300 

51 



2015-02-14 

01:44 

300 

46 


* 60s exposures in dithering mode at 5 positions 


2 OBSERVATIONS 

2.1 The XMM-Newton data 

XSS J1227 was observed twice with XMM-Newton after the 
transition to the radio pulsar state whi ch occurr ed be¬ 
tween end 2012 and beginning 2013 dBassa et al.l 12014) 1. 
The first observation, taken on Dec. 2013 December 29 
(OBSID: 0727961401, hereafter Obs. 1), was analysed by 
IBogdanov et al.l d2014l) . The second pointing was performed 
on 2014 June 27 (OBSID:0729560801, hereafter Obs. 2) (see 
also TablcUl . The former pointing used the EPIC-pn cam¬ 
era JStrrider et al.ll200lh operated in Full Window Mode and 
with a thin optical blocking filter. In such a configuration 
the instrument retains a field of view of 25.6 x 26.2 arcmin, 
and has a temporal resolution of 73.4 ms. During Obs. 2 the 
EPIC-pn was instead set to the Timing fast mode to achieve 
a temporal resolution of 0.03 ms with the aim at detecting X- 
ray pulses at the 1.69 ms rotational period of the NS. In this 
observing mode 2-D imaging capabilities are lost. In both 
observations a Full Window observing mode with a thin fil¬ 
ter w as used for the two EPIC-MOS cameras dTurner et al.1 
l200l|j . giving a time resolution of 2.6 s. Data were filtered 
to remove soft proton flares reducing the effective exposure 
to 34.6 and 4 1.8 ks in Obs. 1 and 2, respectively. As dis¬ 
cussed in IPapitto et al.l d2015l) . the source is not detected in 
the EPIC-pn data obtained in Timing mode during Obs. 2. 
Therefore we focus only on the data acquired with the two 
EPIC MOS cameras and draw comparison between the two 
observations. 

We analyzed data using the XMM-Newton Science Analy¬ 
sis Software (SAS) v.141'0. EPIC MOS source photons were 
extracted from a 40 arcsec wide circle centered at the source 


1 http://xmm.esac.esa.int/sas/ 
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position, to encircle ~ 85% of the source photons. Back¬ 
ground was accumulated from a 50 arcsec wide circular re¬ 
gion far from the source. X-ray photons were reported to 
the Solar system bary centre using th e op tical position of 
the source reported by dMasetti et al.il200q) . XSS J1227 was 
detected by each MOS camera at an average net level of 
0.025(8) cts s _1 and of 0.040(9) cts s _1 in Obs. 1 and 2, re¬ 
spectively. _ 

In Obs. 1 and 2, the Optical Monitor (OM) dMason et alj 
120011 1 was operated in Fast mode, giving a temporal res¬ 
olution of 0.5s, with the U filter (3500-4800 A). Ten ex¬ 
posures were acquired in both observations for a total of 
32.7ks (Obs. 1) and 42.5ks (Obs. 2). Th e OM-U band data 
of Obs. 1 were presented in lBassa et aid (|2014l ). For a com¬ 
parative analysis of the two observations, the two data sets 
were reprocessed and photometry was extracted using the 
SAS task omchain. In a few exposures during Obs. 2 the 
source was not detected and photometry was extracted per¬ 
forming manual detection of the target using the omsource 
task. The source was at an average magnitude of U=19.43(3) 
and 19.92(4) in Obs. 1 and 2, respectively. Correction to the 
Solar system barycentre was also applied. 


2.2 The REM photometric data 

XSS J1227 was observed from 2015 Jan. 20 to Jan. 23 (here¬ 
after Run 1) and from 2015 Feb. 12 to Feb. 14 (hereafter 
Run 2 ) with the 0.6m INAF REM telescope in La Silla, 
Chile dZerbi et al.ll2004l ). The telescope is equipped with the 
ROSS2 earner^ that performs simultaneous exposures in 
the Sloan filters g’, r’ , i’ and z’ and with the REMIR camera 
dConconi et al.l 120041 1 covering simultaneously the near-IR 
in one band. Integration times were 300 s for all optical fil¬ 
ters and a dithering of 5 exposures of 60 s was used for the 
REMIR J-band exposures. The night of Jan. 22 was not pho¬ 
tometric and thus the data are not included here. The total 
coverages were 3.2 h, 3.51 h and 4.1 h in Run 1 and 4.6 h 5.3 h 
5.5 h in Run 2. The log of the photometric observations is 
also reported in Tablc[l] 

The photometric data sets were reduced using standard rou¬ 
tines of IRAF to perform bias and flat-field corrections. Due 
to the low response of the z’ filter the corresponding images 
have not been analysed. For the REMIR observations the 
five dithered images were merged into a single frame. 

For each data set, aperture photometry was performed opti¬ 
mizing aperture radius and sky subtraction was done us¬ 
ing annuli of different sizes. Comparison stars were used 
to check and to correct for variable sky conditions. The 
REM/ ROSS2 and REMIR photometry was calibrated us¬ 
ing the Sloan standard SA 94 242 observed each night whose 
near-IR magnitudes are also tabulated in the 2MASS cat- 
alogu(0. XSSJ1227 was found at g’=18.62(2), r’=18.12(4), 
i’=17.96(2) in Run 1 and at g’=18.75(3), r’=18.28(3) and 
i’=17.91(4) in Run 2. Because the source was barely de¬ 
tected in the individual J band images, image coaddition 
for each night was performed to provide a mean magnitude. 
XSS J1227 was at J=16.98(4) and J=16.92(4) in Run 1 and 


2 http://www.rem.inaf.it 

3 http://www.ipac.caltech.edu/2mass/ 



Figure 1. X-ray orbital modulation observed during Obs. 1 (red 
crosses) and Obs. 2 (blue circles), evaluated in 20 phase bins by 
folding the background subtracted summed light curves observed 
by the two MOS cameras in the 0.3-10 keV band, at the 6.91 h 
orbita l period accor ding to the radio pulsar ephemeris given in 
IPapitto et al.l ll2015l) . Phase 0 corresponds to the passage of the 
NS at the ascending node of the orbit. 
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Figure 2. Orbital modulation of the hardness ratio evaluated as 
the ratio of fluxes in the 2-10 keV and 0.3-2 keV bands in 8 phase 
bins (see Fig.[T}. 

2, respectively. The g’, r’ and i’ light curves were also cor¬ 
rected to the Solar system barycentre. 

3 RESULTS 

3.1 The X-ray orbital variability 

The background subtracted summed MOS light curves in 
the 0.3-10 keV range during Obs. 1 and 2 were folded at 
the or bital pe r iod u sing the radio pulsar ephemeris given by 
IPapitto et alj (|2015l ) where phase 0 is defined as the epoch at 
which the NS passes at the ascending node of the orbit. The 
uncertainties on the determination of the orbital phase with 
these ephemeris is < 1% for the two epochs. The orbital 
modulation observed in Obs. 1 and 2 binned in 20 phase 
intervals, each lasting ~ 1.25 ks is plotted in Fig.|T] 

A minimum count-rate of 0.04 ctss -1 is observed in both 
observations at orbital phase <j> = 0.25, i.e. when the donor 
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star is between the pulsar and the observer. The intensity 
then increases when the pulsar approaches the closest dis¬ 
tance to the observer (</> = 0.75). However, during Obs. 2 
the increase of the flux is steeper and the flux at maximum 
roughly doubles that during Obs. 1. The modulation ampli¬ 
tude, defined as (F max —F m i n )/(F max ±F mill ), evaluated with 
a simple sinusoidal function is 0.25(4) and 0.69(6). Also, a 
dip is observed during Obs. 2 centred at < f> = 0.75 where 
the count rate decreases by a factor ~2. This dip is instead 
marginally visible in Obs. 1. 

An orbital variability is also displayed by the hardness ratio 
evaluated as the ratio of the flux observed in the 2-10 keV 
band to the flux in the 0.3-2 keV band (see Fig.0. The 
orbital modulation was sampled using 8 phase bins in order 
to increase the counting statistics of each point. In both 
observations, the hardness ratio increases by a factor ~ 1.5 
at maximum of the modulation, implying that the spectrum 
gets harder when the pulsar is at inferior conjunction. The 
significance of the orbital variability of the hardness ratio 
is higher in Obs. 2 than in Obs. 1. Modelling the orbital 
profile of the hardness ratio (in Fig. 0 with a constant gives 
Xr = 1-76 and Xr = 2.8 for 7 d.o.f. in Obs. 1 and Obs. 2, 
respectively. This gives a probability that the profiles are 
compatible with a constant of 0.091 and 0.005, respectively. 
The significance of a detection of the orbital variability in 
the hardness ratio is therefore 1.7 and 2.8a in Obs. 1 and 
Obs. 2. In both cases, a sinusoidal fit to the orbital profile 
yields an improvement of the modelling with respect to a 
constant, which has a probability of ~ 5% of being due to 
statistical fluctuations (F-test). We also note that the source 
is significantly detected in all orbital phase bins in the two 
energy bands at a significance larger than 10a. 

3.2 The X-ray spectrum 

For each epoch, the spectra from the two MOS cameras in 
the range 0.3-10 keV were fitted together using the same 
spectral parameters. During Obs. 2 the spectrum is well de¬ 
scribed (xl — 1-18 for 123 d.o.f.) by a power law with index 
r = 1.07 ± 0.08 (all uncertainties on the spectral parame¬ 
ters are quoted at the 90% confidence level). Only an upper 
limit of Nh < 5 x 10 20 cm -2 (90% confidence level) could be 
set on the absorption column. The 0.3-10 keV average unab¬ 
sorbed flux is 7.1±0.4x 1CU 13 erg cm -2 s -1 . The addition of 
a thermal component, either a blackbody (bbodyrad model 
in XSpe c) or a more physical NS atmosphere (nsA model in 
XSpec, (IZavIinll2006i) ). does not give an improvement of the 
modelling. Fixing in NSA the NS mass and radius to 1.4 Mg 
and 10 km, and using the low field option B=0 G, we ob¬ 
tain an upper limit on the 0.3-10 keV flux of such a thermal 
component of 2 x 10 -14 erg cm -2 s _1 , i.e. ~ 3% of the total 
flux. 

We also modelled the spectra during Obs. 2 at orbital phases 
close to the NS superior conjunction (< i> ranging between 
0.1 and 0.4), and to the NS inferior conjunction (</> ranging 
between 0.5 and 1.0). As indicated by the hardness ratio, 
close to the orbital minimum the spectrum gets softer and 
can be described by a power law with index T = 1.304 q’ 2 | 
(Xr = 1-03 for 21 d.o.f.). The 0.3-10 keV unabsorbed flux is 
2.8 ± 0.4 x 10“ 13 erg cm -2 s _1 . A thermal NS atmosphere 
component is not statistically required to model the spec¬ 
trum, and would contribute to at most to 40% of the flux in 



1 10 
Energy (keV) 


Figure 3. The phase-resolved MOS spectra during Obs. 2 (top 
panel) at orbital phases close to the NS inferior conjunction 
(</> = 0.5 — 1.0, blue points) and close to the NS superior con¬ 
junction (</> = 0.1 — 0.4, red points). Middle and bottom panels 
show residuals with the best fit absorbed power law model, shown 
as a black solid line. 

the same energy range (l.ljjg'® x 10~ 13 erg cm -2 s _1 ). Close 
to the NS inferior conjunction the spectrum is described by 
a harder (T = 1.0 ± 0.1) power law (y 2 = 0.93 for 73 d.o.f.) 
that emits a 0.3-10 keV flux of 1.08 ±0.08 x 10 -12 erg cm -2 
s _1 . In this case, a thermal component would contribute to 
at most 13% of the total flux. The fits to the MOS spec¬ 
tra and those at the two phases using the simple absorbed 
power law are shown in Fig. [3] 

The fit to the MOS s pectra during Obs. 1 is consistent with 
the results found bv lBogdanov etakl J20141 ) for the EPIC- 
pn spectrum, the power law index being F = 1.2 ± 0.1 
(x% = 0.87 for 79 d.o.f.) and an upper limit to the hy¬ 
drogen column density Nh < 6.8 x 10 2 ° cm -2 (90% confi¬ 
dence level). The unabsorbed flux in the 0.3-10 keV range 
is 4.4 ±0.1 x 10“ 13 ergcnU 2 s . Due to the lower sensitiv¬ 
ity of the MOS cameras with respect to the EPIC-pn, we 
are unable to conclude on the detection of the NS thermal 
emission. Here we note that the spectral fits to the EPIC-pn 
data d o not i mp rove but not exclude a thermal component 
llBoedanov et al.ll2014l ). 

Comparing the power law indexes during the previous disc- 
accretion state (dM10,dM13) (r = 1.7) the spectra are 
harder during the disc-free states in Obs. 1 and 2 and the 
absorption contribution is consistent with that previously 
determined in dM10 and dM13. The luminosity in the 0.3- 
lOkeV range, assuming _d=l.4 kpc as derived from the ra¬ 
dio DM measures dRov et al] l2015| ). results to be 1.0-1.7 
xl0 32 ergs _1 at the two epochs, one order of magnitude 
lower than in the disc-accretion state. 

3.3 The UV/optical orbital modulation 

We folded the UV photometric data of Obs. 1 and 2 at the 
orbital period using the same ephemeris adopted for the X- 
ray band. The orbital modulation in the UV band at the 
two epochs is reported in Fig. [4] The UV flux is minimum 
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Figure 4. Orbital modulation observed in the OM U band during 
Obs. 1 (red crosses) and Obs. 2 (blue circles), evaluated in 20 
phase bins. 

at 0 ~ 0.25 when the source shows roughly the same mag¬ 
nitude (U~20.4) in both observations. The UV flux reaches 
a maximum at 0 ~ 0.75 where it increases by 1.43(4) mag 
and by 0.72(9) mag in Obs. 1 and 2, respectively. The slight 
decrease at 0 ~ 0.9 in Obs. 2 is not significant from inspec¬ 
tion of the unfolded light curve. The UV and X-ray band 
modulations are thus in phase, but the amplitudes in the 
two bands behave in opposite way. In Obs. 1 the X-ray am¬ 
plitudes are larger than in Obs. 2, but the UV amplitudes 
are lower in Obs. 1 than in in Obs. 2 (see Fig.[l]). 

We folded the optical g’, r’ and i’ REM photometry of 2015 
Jan. (Run 1) and Feb. (Run 2) at the orbital period using 
the same ephemeris. In Fig. [5] we report the orbital mod¬ 
ulation observed at the two epochs together with the g’-i’ 
colour index. Likewise the X-ray and U bands, the modula¬ 
tions in the optical/near-IR bands show a minimum centred 
at 0=0.25 and a maximum at 0=0.75. As evaluated with si¬ 
nusoidal fits the modulation amplitude is Ag'=0.84(2) mag, 
Ar'=0.74(l) mag and Ai'=0.56(l) mag in Run 1, and it is 
Ag'=0.72(2) mag, Ar'=0.52(l) mag and Ai'=0.40(l) mag in 
Run 2. The source is also on average brighter in g’ and r’ 
by 0.1 mag in Run 1 than in Run 2. A comparison with the 
optical level observed in 2013 Dec. (dM14), the source was 
at about the same g’ magnitude as in Run 2, with a modula¬ 
tion of similar amplitude. The modulation is also colour de¬ 
pendent and exemplified by the g’-i’ index (bottom panels of 
Fig-EJ, indicating a reddening at orbital minimum when the 
donor star is at inferior conjunction. The average values are 
g’-i’=0.66(6) and g’-i’=0.84(4) in Run 1 and 2, respectively. 
The source thus appears redder when on average fainter. 

The large amplitude modulation and strong colour variabil¬ 
ity confirm the previous finding that the donor star suf¬ 
fers irradiation (dM14). Similarly, we applied the Night¬ 
fall cod(Q to fit simultaneously the g’, r’ and i’ light curves 
with the aid of the donor star radial velocity curve presented 
in dM14. We adopted a blackbody temperature of 5000 K 
as starting point for the unheated star and fixed the max¬ 
imum allowed temperature (500 000 K) for the irradiating 
star treated as a point source. The donor star is assumed to 
fill its Roche lobe which is consistent with the occurrence 

4 The Nightfall code is available at 
http: //www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall. 


of radio eclipses for ~ 40% of the orbit llRov et al.ll2015l ). 
We left the mass ratio q, the binary inclination i and the 
donor star temperature T 2 free to vary. We derive for Run 1 
q = 0.09 ± 0.08, i = 53^° and T 2 =6000±500 K heated up 
to 6500 K (xr=9.5 for 362 d.o.f.) (uncertainties are at lcr 
confidence level) and for Run 2 q = 0.191q jg, i = 50l^g° 
and T2=5900±400K heated up to 6300 K (Xr=15.3 for 482 
d.o.f.). The large \r are due to the small errorbars of the 
photometric measures. The g’, r’ and i’ light curves show 
asymmetries in the portion ri sing to t h e ma ximum but the 
inclusion of cold or hot spots IIRomani et alj|2012i ) does not 
improve the fits. Similar behavio ur is als o found in_the black 
widow binary PSR J1311-3438 bv lRomani et all (120151 ). sug¬ 
gesting that the donor face is not fully heated by a point 
source and that the asymmetric shock may provide at least 
some of the heating. We also used the U band photometry of 
Obs. 1 together with the donor star radial velocities and find 
q = 0.16j)Q Qg, i = 58^2° an d T 2 =5500±300 K heated up to 
60001< (x^=6.2 for 88 d.o.f.). The addition of the g’ and r’ 
photometry acquired in dec. 2013 (dM14) does not improve 
the results because of the large scatter and the sparse orbital 
coverage of the data. In all three data sets the mass ratio q 
is c onsistent w ithin errors with that derived from radio data 
by IRqv et ahi (|2015i ). q=0. 194(3). We note that the large 
filling factor of the donor implies an aspherical shape for 
the star that causes strong temperature (brightness) varia¬ 
tions over the surface and the model code provides a mean 
temperature of the heated face of the secondary. The g’, r’, 
i’ fluxes at the peak of the modulation observed in Run 1 
and 2, when dereddened for intervening absorption Eb-v = 
0.11 (dM10), are consistent with a blackbody with temper¬ 
atures similar to those derived from the light curve models, 
6600±200 K and 6300±300 K, respectively. The blackbody 
normalisation decreases by a factor of 1.17±0.01 between 
Run 1 and 2. We therefore confirm the findings in dM14, 
that the donor star temperature is compatible with that 
of a G5 star seen at inferior conjunction and it reflects a 
F5 spectral type at superior conjunction. Since the binary 
parameters q and i cannot explain the variability in the am¬ 
plitudes observed at the different epochs, we also conclude 
that the variations are due to changes in the visible area of 
the irradiated face of the donor star (see also Sect. 4). In the 
following we will assume 0.11 < q < 0.26 and 46° < i < 65°. 
For a NS mass of 1.4 Mq, the mass ratio gives a mass of 
donor star 0.15 < M 2 < 0.36. 


4 DISCUSSION 

We here discuss the X-ray and UV/optical orbital variability 
of XSSJ 1227 during the current disc-free state. 

4.1 The X-ray and optical variabilities 

The X-ray flux is strongly modulated and it is maximum 
close to the inferior conjunction of the NS. The mini mum 
occur s at the same phases of the radio eclipses (iRov et alJ 
120151) . which are produced by the orbital motion of the sec¬ 
ondary star. X- ray orbit al modulations are also observed in 
other redbacks dRoberts et al.ll2015l ). The X-ray spectrum as 
well as those at minimum and maximum of the modulation 
ilare non-thermal. A dominant non-thermal spectrum may 
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Figure 5. From top to bottom: The orbital modulation in the optical g’, r’ and i’ bands and the g’-i’ index as observed in Run 1 (left) 
and Run 2 (right). Single nights are reported with points of different colours and the average of the three nights as black points evaluated 
in 30 phase bins. The black lines report the best fit Nightfall curves (see text for details). 


suggest that the X-ray emission originates in an intrabinary 



the orbital cycle, suggests that the X-ray emission is ex¬ 
tended and originates close to the donor star and presum¬ 
ably at the LI point, supporting an origin in the intrabi¬ 
nary shock. Si milar interpretation is given for othe r eclips¬ 
ing redbacks (iBogdanov et al.l l201ll ; ITT et al.ll2014l ). where 
the X-ray modulation is caused by the eclipse of the in¬ 
trabinary shock emission by the donor. The eclipse depth 
and duration depend on the ratio between the donor star 
radius and the intrabinary shock separation and on the bi¬ 
nary inclination. The shock location or apex is defined as 
r s /a ~ 1/(1 + r? 1 / 2 ) (I Arons fe Tavanil 11993lf . where a is the 
binary separation and r\ ~ Mv w c/L s d is the momentum 
ratio between the stellar wind and the pulsar wind. Unfor¬ 
tunately, this quantity cannot be determined because the 
mass loss rate during the radio pulsar state is unknown 
and thus also the shock apex. For what follows we assume 
r s ~ Rl.ns ~ 1.1 — 1.2 x 10 11 cm adopting the mass of the 
components in Sect. 3.3. For an isotropic pulsar wind the ge¬ 
ometric fraction of the wind stopped by the shock, assumed 
to be of the same size of the Roche-lobe filling donor star, 
is f s ~ 0.5 (1 — cos ft ) ~ 0.02, where fl = atari (R2/Rl,ns) 
and R 2 = Rl2 = 0.462 (q/1 + q) 1/3 a ~ 3.2 - 4.2 x 10 10 cm 
for the range of q found in Sect. 3.3. This fraction is sim ilar 
to those inferred in other redbacks llRoberts et al.ll2015l ). 


The shock is expected to emit synchrotron radiation. To 
esti mate the synchro t ron lu mino sity from the shock we fo l¬ 
low lArons fe Tavanil (1 19931 ) and lKennel_jg L ^lqronitj| ll 198^1. 
Similar arguments were also used by Stappers et al.l fl2003l ) 
for the black widow binary PSR B1957+20 and first in¬ 
voked for the redback bi nary PSR J0024-7 2 04W in the glob¬ 


ular cluster 47 Tuc by IBogdano v et al. l (120051) and later 
also for PSR J1023+0038 bv IBogdanov et al. (1201 ll) . The 


X-ray luminosity emitted by the shock depends on both the 
post-shock magnetic field strength and the ratio of mag¬ 
netic and particle energy density a. This fraction is not 
known but could be as low as cr=0.003 in winds domi¬ 
nated by kinetic energy, like the Crab pul sar, or as high 
as a » 1 for strongly magnetised winds dStappers et al.l 
120031) . Given the small separation between the pulsar and 
the companion, the intrabinary sh ock is expected to be 
located in a strong magnetic field dArons fe Tavanil Il993l : 
IStappers et al.l 120031 ). The magnetic field upstream of the 
shock is Bi = [(cr/(l + ct) (E/c f p r 2 )] 1 ^ 2 , where f p is the pul¬ 
sar isotropic factor and E is the spin-down energy of the 
pulsar. In the two extreme c ases, assuming a pulsar energy 
loss of E ~ 9 x 10 34 ergs _1 dRov et aid 120151) . an isotropic 
wind f p =l and the shock at r s ~ Rl.ns ~ 1.2 x 10 11 cm, we 
obtain Bi ~ 0.8 G and ~ 16 G, respectively. The post-shock 
magnetic field is B 2 ~ 3Bi ~ 2.5 G for a « 0.01 and it 
is B 2 ~ Bi ~ 16 G for high- <7 values dKennel fe Coronitil 
119841) . For simplicity we assume that the observed X-ray 
luminosity totally originates via synchrotron emission from 
the shock. This requires that the electrons are relativistic 
with 7 ~ 2.5 x 10 5 (e/B 2 ) 1 ' 2 , where e is the observed photon 
energy, ~ 10 keV. Thus 7 ~ 2 — 5 x 10 5 in the case of mag¬ 
netic and kinetic energy dominated winds, respectively. The 
expected synchrotron luminosity is L sync ~ f syn c f 7 f s L s< j, 
where L s d is the spin-down luminosity, f 7 is the fraction 
of spin down energy flux that accelerates electrons with 
Lorentz factor 7 at lOkeV, and f sync = (1 + t syn c/tfl 0 w ) _1 
is the radiative efficiency of the synchrotron emission. Here 
ts yn c is the radiative synchroton loss time and tn ow is the 
time of residence of particles radiating in the shock. For the 
assumed shock size of the order of the donor star diameter 
tflow ~ (3r/c) ~ 6- 9 s for a « 1 and tfl ow ~ (r/c) ~ 2- 
3 s for high-a flows dKennel fe Coronitilll984l ). The radiative 
loss time t sync = 5 x 10 8 /( 7 B 2 ) is ~ 10 s and ~160s for 
a >> 1 and a << 1, respectively. Thus f sync ~ 0.04 — 0.05 
for kinetic dominated winds and f sync ~ 0.17-0.23 for mag- 
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netic dominated winds. Using the geometric factor f a ~ 0.02 
derived above and the observed 0.3-10 keV X-ray luminos¬ 
ity at orbital maximum of 2.5 x 10 32 ergs -1 , we obtain 
f 7 ~ 0.60 — 0.81 for the highly magnetic case and unphysical 
values f 7 ~ 3 for the kinetic dominated wind case. These 
large values reflect the geometrical uncertainties of the pul¬ 
sar wind f p and the fraction stopped by the shock f s , indi¬ 
cating that the pulsar wind is not isotropic and hence the 
f a > 0.02 should be regarded as a lower limit. The unphysi¬ 
cal values for f 7 in the low-cr limit also imply that the wind 
must be magnetically domin ated. Similar conclusi ons were 
drawn for PSR J1023+0038 llBogdanov et al .11 201 lh . 

The X-ray spectrum hardens at orbital maximum with 
T ~ 1.0, whilst at minimum T ~ 1.3. A hardening at orbital 
maximum is a new observational aspect in r edbacks, as it has 
been only detected in PSR J2129-0429 llHui et al .1120151 ) so 
far. It is instead observed in other more massive gamma-ray 
binaries such as LS 5039 (Li et al.|[201lh . The hardening can¬ 
not be due to photoelectric absorption since it would occur 
at orbital minimum. It could suggest the presence of an addi¬ 
tional contribution to the X-ray emission. In spin-powered 
X-ray pulsars, a fraction ~ 10” 4 of the spin down lumi- 
no sity is converted in X-rays. Using the empirical relation 
bv iPossenti et al.1 ll2002l b for L s d ~ 9 x 10 34 ergs 1 an X- 
ray luminosity of magnetospheric origin ~ 3.5 x 10 31 ergs” 1 
is expected in the 2-10keV range. This is not much dif¬ 
ferent from that observed at minimum in the same energy 
range: Lx,min ~ 4.9 x 10 31 ergs” 4 . This emission should 
be pulsed at the spin period of the NS. Unfortunately the 
source is n ot detected in t he EP IC-pn timing data during 
Obs. 2, but IPapitto et al.l (|2015l l derive an upper limit on 
the 0.5-10 keV pulsed luminosity of ~ 1.6 x 10 31 ergs _1 . 
The limit on the NS thermal emission at minimum is esti¬ 
mated as ~ 0.4Lx,min ~ 2.5 x 10 31 ergs” 1 in the 0.3-10keV 
range (sect. 3.2). Furthermore, at X-ray orbital minimum 
the shock could still contribute, because of its large size and 
the energy distribution of the electron population may be 
softer in the outer shock regions. Therefore, it is not pos¬ 
sible to conclude whether the observed softening at orbital 
minimum is due to the extended shock, to the NS or both. 
Furthermore while only a hint of a dip superposed on the 
broad maximum could be detected in Obs.l, this feature is 
observed when source is brighter in Obs.2 (see Fig.|T]). A 
doubl e-peaked maximum is is also prese nt in PSRJ2129- 
0429 llHui et al.ll2015l : iRoberts et al.ll2015l l at similar orbital 
phases. In PSR J1023+0038 an enhancement of e mission 
is observed at eclipse egress llBogdanov et al.l 1201 il l. These 
variations are interpreted as due to Doppler boosting due 
to the finite velocity of the shock wind in the direction of 


thej 

aost-shock flow 1 

Arons & Tavani 1993|; Staooers et al. 

2003 

: Bogdanov et al 

201l!: Huang et all 2013). If the shock 


wraps the donor star, the plasma accelerates as it flows 
around the companion and, around inferior conjunction of 
the donor star, the electrons are accelerated primarly in the 
direction of the observer. At this phases an increase of flux 
would be expected, but due to obscuration by the donor star, 
the enhancements would occur before and after the eclipse. 
For PSR J1023+0038 the absence of enhancement prior the 
eclipse ingress was interpreted as due t o a weaker out¬ 
flow toward the leading edge of the eclipse llBogdanov et al.l 
201l|). For post- s hock velocities in the range c/3 — c/3 1 ^ 2 , 
Arons fc Tavanil (Il993l l predict brightness variations by a 


factor of ~1.3-2.2, respectively. The flux variations at the 
dip of a factor ~ 2 observed in Obs. 2 is consistent with this 
prediction. 

In Obs. 2 the X-ray peak flux at maximum is ~2 times 
larger than on Obs. 1. If the mass loss rate and thus the 
luminosity increases, then the momentum ratio rj increases 
and the shock apex moves inward the NS Roche lobe. How¬ 
ever, a change by a factor of ~2 in the mass loss rate will 
marginally affect the apex radius, but it may produce an 
increase in size of the shock. This would imply an increase 
of the time of residence of electrons tfj ow in the shock and 
of the fraction of the pulsar wind intercepted by the shock 
(f s ). It could be also possible that the relativistic plasma 
accelerating as it flows around the companion may reach 
larger velocities giving rise a larger Doppler boosting effect 
as observed in Obs. 2. 

The orbital modulation in the UV/optical is strong, amount¬ 
ing to ~ 0.7-0.8 mag in the optical and up to ~1.4mag in 
the UV. The optical colours along the orbital cycle indicate 
a reddening at inferior conjunction of the donor star. The 
light curves are well describe a G5 type donor star filling 
its Roche lobe and heated up to ~ 6500 K, although the 
modelling considers a point-like irradiating source at the 
NS and not at LI. Thus the donor star is heated by high 
energy radiation from either the shock and/or by the pul¬ 
sar wind. We here recall that irradiation of the secondary 
was also found to occur before the state transition (dM14). 
As also noted in dM14, the donor is not expected to sub¬ 
stantially cool, because the thermal relaxation time scale 
is of the order of the Kelvin-Helmhotz time of the convec¬ 
tive envelope, tkh ~ (M 2 /R 2 L 2 ). The observed changes in 
the amplitudes of the UV/optical orbital modulation would 
then be due to changes of the visible area of the heated 
face of the companion rather than to temperature varia¬ 
tions. Thus, the opposite behaviour seen in the X-ray and 
UV modulation amplitudes could be explained in terms of a 
larger size of the shock, producing a stronger X-ray mod¬ 
ulation, and at the same time preventing the irradiated 
face to be totally visible, yielding to a weaker modulation 
in the UV and optical bands. The heating effect on the 
donor star, Li rr = 4 7r R| a Tf rr ~ 1.2 — 2.2 x 10 33 ergs” 1 
for Ti rr ~ 6500 K0, is larger than that estimated for the 
shock ~ 3x 10 32 ergs” 1 in the 0.3-10 keV range. However, if 
the X-ray spectrum extends well beyond 10 keV, as observed 
by NuSTAR in PSR J1023+0038 with no spectral break up 
to 79keV (iTendulkar et al]|2014lb the X-ray luminosity at 
orbital maximum in the 0.3-100 keV range can be about one 
order of magnitude higher than that measured in the XMM- 
Newton EPIC range. In this case the shock could be totally 
responsible for the heating of the donor star. If instead we 
consider the pulsar spin-down power, we infer that ~ 0.45 of 
the spin-down luminosity goes into the heating of the donor 
star. Observations at higher energies will clarify this issue. 


5 A detailed and rigorous treatmen t of the irradiation of the 
donor star can be found in lRitter et ali (12000 1. The effect of ir¬ 
radiation is to suppress the temperature gradient in the outer 
stellar envelope and blocking the outward transport of energy 
through the face. Thus the emergent flux from the irradiated face 
is: F; rr = r7T/. r = WF; nc , with W~0.5, the albedo for a star in 
convective equilibrium. 
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Figure 6. The broad-band spectral energy distributions (SED) of XSS J1227 during the rotation-powered state (left panel) and during 
the sub-luminous pre-transition state (right panel). The post-transition SED is constructed with the X-ray spectra of the two MOS 
cameras (blue and magenta), the UV O M and optical REM photometry (red squares), the radio measures reported by |Rove^rlJjj2015|) 
(blue crosses) and the upper limit bv lBassa et al.l ((2014 ) (blu arrows) and the Fermi- LAT (light blue) fluxes reported in IXin^fcWang 
1120141 ') together with the corresponding fitted model (see text for details). The pre-transit i on SE D is taken from dM13 but substituting 
the Fermi -LAT fluxes reported in the second release catalogue with those in lXing fc Wanei ( 2014 ) together with the corresponding model 
fit. 


4.2 The spectral energy distribution 

We constructed a broad-band spectral energy distribu¬ 
tion (SED) from radio to gamma-rays during the rotation- 
powered state, shown in the left panel of Fig QJl We use d the 
radio flux at 607 MHz as measured by I Roy et al.l d2015i ) and 
the upper limits ob tained at 1.4, 5.5 and 9 GHz reported by 
iBassa et al.l (j2014l ). These are broadly consistent with the 
spectral slope expected fr om a radio pulsar with a power 
law index -1.7 dRov et al.il2015l ). For the UV/optical range 
we use the average U and g’, r’, i’ and J band magnitudes 
obtained in 2014 and 2015 respectively, corrected for inter¬ 
vening absorption Eb-v = 0.11 (dM10) and overlay the best 
fit blackbody model with Tbb = 5200±300 K. The unab¬ 
sorbed average X-ray spectra in the 0.3-10 keV range of the 
two MOS cameras during Obs 2 are reported together with 
the best fit power law with T=1.07. We also add the Fermi- 
LAT s pectrum after the state transition from lxing fe Waned 
d2014l) . together with a power law with an exponential cut¬ 
off with index T=1.8 and energy cutoff E c = 2 GeV obtain ed 
by the same authors. Here we note that also I Johnson et all 
d2015l) analyse the Fermi -LAT data obtaining comparable 
values. In the right panel of Fig[S]we also report the SED ob¬ 
tained during the previous LMXB state presented in dM13, 
substituting the Fermi -LAT fluxes, obtained at that time 
from the second r elease of the source ca talogue (2FGL), with 
those reported bv IXing fe Wan j d2014l ) during the same pre¬ 
transition state together with their best fit power law expo¬ 
nential cut-off with index F=2.13 and E c = 6 GeV. Also 
in this case the fit parameters are compatible withi n un¬ 
certai nties with those recently obtained by Ijohnson et alJ 
ll2015l) analysing the same period. 

The SEDs at the two epochs substantially differ from each 
other implying very different contributions. I 11 the post¬ 
transition state, the radio flux is a factor ^6 lower than 
that observed in the disc state and the slopejs consis¬ 
tent with that of a radio pulsar (|Rov et ali l2014ll . The 


UV/optical/near-IR emission totally originates from the 
donor star, whilst an additional contribution from the disc 
is required in the LMXB state (dM14). The X-ray emission 
in the radio pulsar state is harder than that before the tran¬ 
sition (r ~ 1.7, dM13). In the LMXB state the X-ray and 
gamma-ray emission s have been successfully modelled by 
a pro peller scenario llPapitto et al.ll2014l ; IPapitto fe Torres 
1201511 with marginal accretion from the NS, while in the 
rotation-powered state the X-ray emission mainly originates 
from the intrabinary shock. The post-transition gamma- 
ray emission appears instead to originate in the pulsar 
outer magnetosphere (Ijohnson et al.l 120151 ) from the ob¬ 
served partially alignment of the Fermi -LAT and radio 
spin light curves and the lack of an orbital modulation at 
GeV energi cfl The gamma-ray luminosity was found by 
Ijohnson et al] (120151) to be ~ 5% of the pulsar spin-down 
energy. Although the lack of X-ray coverage extending above 
10 keV cannot allow us to draw conclusions, the present data 
suggest that the X-rays and high energy gamma-rays are not 
linked. 


5 CONCLUSIONS 

We have analysed the X-ray and optical emissions of 
XSSJ1227 during the ongoing rotation-powered state into 
which the source entered at the end of 2012/beginning 2013. 
We here summarize the main results: 

• The X-ray luminosity in the 0.3-10 keV range during 
the radio pulsar state is ~ 1.0 — 1.7 x 10 32 ergs^ 1 , one or¬ 
der of magnitude lower than that in the previous accretion- 

6 The claim of a weak orbital modulation by Ixing fc Wand 
j20lj) during the r otation-powered state is not confirmed by 
Ijohnson et all ll2015l ) 
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powered state. It is also found to vary by a factor ~1.6 
between two observations taken 6 months apart. 

• We confirm the str ong orbital X-ray and UV and optical 
modu lations found by iBogdanov et al.l (l2014h , iBassa et al.l 
112014) and dM14, respectively. 

• The X-ray modulation amplitude is found to vary by a 
factor of ~2.7 betweem the two epochs. The variability is 
consistent with X-rays originating at the intrabinary shock 
between the pulsar and donor star, which is eclipsed due to 
the companion orbital motion. Doppler boosting could be 
important to shape the X-ray modulation when the source 
brightens. 

• The X-ray spectrum is non thermal (T=1.07) and it 
is harder than that observed during the accretion-powered 
state and consistent with synchrotron emission from the in¬ 
trabinary shock. The spectrum also softens at orbital min¬ 
imum, a new feature seen only in one other redback so far. 
This could be due to a softer contribution of the outer re¬ 
gions of the shock and/or of the NS. 

• The UV/optical/near-IR emission is modulated at the 
orbital period due to irradiation of the companion star. We 
infer a 5500 K temperature for the unheated hemisphere of 
the donor and a ~6500 K temperature for the irradiated face. 
We also derive a binary inclination 46° < i < 65° and a mass 
ratio 0.11 < q < 0.26. For Mns = 1.4Mg, the donor mass 
is 0.15 < M 2 < 0.36 Mq, confirming the redback nature of 
this binary. 

• The UV/optical orbital modulation also changes in am¬ 
plitude with time. The simultaneous UV and X-ray modu¬ 
lations behave opposite, the UV flux and amplitude being 
larger when the X-ray flux and amplitude are lower. This be¬ 
haviour could be due to small changes in the mass flow rate 
from the donor star which affect the size of the intrabinary 
shock. 

• The broad-band SED from radio to gamma-rays in the 
disc-free state is different from that observed during the 
LMXB state and composed by multiple contributions. The 
radio and GeV emissions are compatible with an origin from 
highly accelerated particles in the pulsar magnetosphere, 
while the X-rays are dominated by the intrabinary shock. 
The UV/optical and near-IR fluxes instead originate in the 
donor star. Whether the companion is heated by the intra¬ 
binary shock or by the pulsar spin-down power needs to be 
assessed with a higher energy coverage in the X-rays. 
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